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ABSTRACT 

We present deep, wide-field imaging of the M51 system using CWRU’s Burrell Schmidt telescope at KPNO 
to study the faint tidal features that constrain its interaction history. Our images trace M51 ’s tidal morphology 
down to a limiting surface brightness of fniim ~30 mag arcsec 2 , and provide accurate colors (<r^ _ y < 0.1) 
down to jig ~ 28. We identify two new tidal streams in the system (the South and Northeast Plumes) with 
surface brightnesses of /ig = 29 and luminosities of ~ 1 0 6 Lq, ; . While the Northeast Plume may be a faint 
outer extension of the tidal “crown” north of NGC 5195 (M51b), the South Plume has no analogue in any 
existing M51 simulation and may represent a distinct tidal stream or disrupted dwarf galaxy. We also trace 
the extremely diffuse Northwest Plume out to a total extent of 20' (43 kpc) from NGC 5194 (M51a), and 
show it to be physically distinct from the overlapping bright tidal streams from M51b. The Northwest Plume’s 
morphology and red color (B-V = 0.8) instead argue that it originated from tidal stripping of M5la’s extreme 
outer disk. Finally, we confirm the strong segregation of gas and stars in the Southeast Tail, and do not detect 
any diffuse stellar component in the H I portion of the tail. Extant simulations of M51 have difficulty matching 
both the wealth of tidal stmcture in the system and the lack of stars in the H I tail, motivating new modeling 
campaigns to study the dynamical evolution of this classic interacting system. 

Subject headings: galaxies:individual(M51), galaxies:interactions, galaxies:evolution, galaxies:spiral 


1. INTRODUCTION 

Galaxy interactions form the basis of the current hierarchi¬ 
cal accretion paradigm for galaxy evolution, and have been 
linked to the formation of bars, spiral structure, and induced 
star formation. However, all these features may also result 
from purely secular processess (see the review by Athanas- 
sou!a1 [20T0| , entangling the various evolutionary effects in 
individual systems. Understanding interaction histories of 
nearby galaxies, from which much of our detailed understand¬ 
ing of physical processes comes, is thus particularly impor¬ 
tant. 

The galaxy pair M51 is one of the best known nearby in¬ 
teracting systems, comprised of the grand-design spiral NGC 
5194 (M51a) and its early-type SB0 companion NGC 5195 
(M51b). Deep imaging in the optical and infrared (e.g. Burk- 
|head|[l978j |Smith et al. |1990[ |Mutchler et al.| 1 2005 i cou- 


pled with kinem atic studies in o ptical (e.g. Tully|1974 
rell et al. 2003 1 and radio (e.g. Appleton et al. [T086 


Dur- 

R’ots 


et al.||1990|>, have revealed extended tidal streams and com¬ 


plex kinematics, including an apparently counter-rotating tail 
of gas on its southern end (jRots et al.|1990 1 . 

M51’s proximity (7.5 Mpc; [Ciardullo et al.||2002| 

^ .ed 


Bose & 

|Kumar||2014]l and well-defined spiral structure has made it a 
frequent target for studies of the interplay between spiral dy- 
namic s, star formation, and the atomic and molecular ISM 
(e.g. |Calzetti et al.|2005[|Hughes et al ,|20 1 3[ ). If much of the 
star formation present in the system is interaction-induced, 
models become an important means of constraining the phys¬ 
ical mechanisms and timescales involved, and hence, for ex¬ 
ample, the calibration of star formation tracers in galaxies 
(e.g. ICalzetti et al. 2005 | and determination of generalized 
star formation laws (e.g. |Kennicutt et al.|2007| l. 

This has made the pair a favorite subje ct of interaction mod- 
els (e.g. Toomre & Toomre 1972 Howard & Byrd|1990[ Salo 
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|& Laurikainen|[2000l |Durrell et al.||2003l |Dobbs et al.||2010] 

and many others), constrained by the morphology and kine¬ 
matics of the observed tidal features. These modeling cam¬ 
paigns have led to two classes of models for the system: one 
involving a single fly-by passage (Toomre & Toomre 1972 


Hernquist|1990 Durrell et al.|2003| l, and a sec ond involving 

multiple passes on a more tightly bound orbit dSalo & Lau- 


tightly 


rikainen 2000; Dobbs et al. 2010). Discriminating between 


these scenarios can constrain both the timescale and strength 
of the perturbation that set up the grand-design spiral and in¬ 
duced star formation in the M51 system. 

Deeper observations can prove useful in locating additional 
interaction signatures to further constrain the models. Partic¬ 
ularly in M51, in which both galaxies have vastly different 
star formation rates ( [Lee et ah] 201 1| |2012[ >, broadband col¬ 
ors may be used to help determine the origin of many inter¬ 
action signatures, and add additional constraints from stellar 
populations. Combining these observations with the kinemat¬ 
ics and spatial distribution of the gas in the system can also 
greatly narrow down the available parameter space for the dy¬ 
namical models. In this Letter , we present deep multiband 
optical imaging of M51, taken as part of a survey of low sur¬ 
face brightness features of nearby galaxies, using the 0.6/0.9m 
CWRU Burrell Schmidt Telescope at Kitt Peak National Ob¬ 
servatory (jWatkins et al.|2014). We report on the morphology 
and colors of both known and newly-discovered tidal features 
which will help constrain the dynamical evolution of the M51 
system. 

2. OBSERVATIONS AND ANALYSIS 

We give a brief account of our observation and data reduc¬ 
tion procedures here; more details can be found in Watkins et 
|ak| ( |2014) > and references therein. We observed M51 in spring 
2010 and spring 2012, with a setup identical to that described 
i Watkins et al. (2014|. We observed only on moonless pho 


tometric nights, in two filters: modified Johnson B (31x1200s 
exposures) and Washington M (39x900s exposures). We gen¬ 
erated dark sky flats using a comparable number of offset 
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Figure 1. Left: a subset of our B band mosaic, rescaled in intensity for different surface brightness ranges: p < 24.6, 24.6 <p/y< 26.5. andp/y> 26.5. At the 
lowest surface brightnesses, the image has been median binned to enhance faint structures. Right: as left, marking regions detailed in Tableftland discussed in 
the text. Each image is 39' (84 kpc) on a side. 


Table 1 

Photometry of Tidal Debris 


Region 

H B,peak 

M b 

B-V 

10 6 Lo, s 

Origin 

Gas Content 

West Stream - North 

25.43 

-14.3 

0.75 

45 

M5 lb 

Gas Poor 

West Stream - South 

25.46 

-14.2 

0.73 

39 

M5 lb 

Gas Poor 

East Stream - North 

25.78 

-14.1 

0.73 

39 

M5 lb 

Gas Poor 

East Stream - South 

25.63 

-13.9 

0.72 

30 

M51b 

Gas Poor 

Southeast Tail 

25.80 

-14.8 

0.64 

68 

M51a 

Gas Rich, offset 

Northwest Plume 

27.59 

-14.5 

0.81 

101 

M51a? 

Patchy? 

Northeast Plume 

29.0±0.3 

-11.2 

— 

4 

M51a? 

Gas Poor 

South Plume 

29.2±0.7 

-11.8 

— 

8 

Unknown 

Gas Poor 


Note. — Absolute photometric errors are cr B ~ 0.04 and (Tb-v ~ 0.05. 


night sky images in each filter (Watkins et al. 20141. Pho¬ 
tometric c onversion to sta ndard Johnson B and V used obser¬ 
vations of |LandoIt| ([1992 ) standar d fields to determi ne color 
terms, and on-frame SDSS (DR8; |Aihara et al.|2011[ > stars to 
derive direct photometric zeropoints. 

Data reduction consisted of standard reduction procedures 
and modeling and subtraction of reflections and extended 
wings of bright stars (see Slater et al7| |20Q9[ ). Following sky 
subtraction, each image was spatially registered and scaled to 
a common photometric zeropoint before median combining 
all images into the final B and V mosaics. We then masked the 
final mosaics of bright discrete sources and spatially rebinned 
and median filtered them in 9x9 pixel blocks to enhance faint, 
diffuse features. Our final B image and our colormap of M51 
are shown in Figures[]]and[2j respectively. Our limiting mag¬ 
nitude is /d B ~ 30. 

In the quantitative photometric analysis that follows, we 
use polygonal apertures to define regions of interest in M51 
(see Figure [T] and Table [TJ. In these regions, contamination 


from background sources can be significant, so we subtract 
off a loc al background flux (see |Rudick et al. 2010; Mi hos et| 
al.|2013 for details). At extremely low surface brightnesses, 


this background contamination dominates the uncertainty, and 
is non-Gaussian. At surface brightnesses of /j b < 27.5, the 
relative uncertainties are small, but at extremely low surface 
brightness (fi B > 28.5) they become significant and preclude 
measurement of meaningful colors. All values have also been 
corrected for the local foreground extinction (A B = 0.126, 
E(B-V) = 0.031;|Schlafly & Finkbeinerj201if. 

In searching for tidal features in deep surface photometry, 
an additional source of contamination is scattered light from 
foreground Galactic dust, which can mimic diffuse structures 
(see |Rudick et al.||20f0f . However, such dust radiates in 
the mid- and far-infrared, and a comparison of our ima ges 
with the WISE 1 2/urn and IRIS lOO/rm images (Meisner & 
Finkbeiner |2014[ |Miville-Deschenes & Lagache||2003] ~ re¬ 


spectively) shows no evidence for strong dust contamination 
cospatial with the diffuse features we identify in Figure [T] (to 







































Extended tidal debris of M51 


3 



Figure 2. B-V colormap of the M51 system. The image spans 36' X 27' (78 kpc X 58 kpc). 

limits of 3.58 MJy sr- l /l.llMJy sr 1 in WISE/IRIS). 

3. RESULTS 

3.1. Morphology and color of tidal features 

An examination of Figure [T] shows a wealth of tidal de¬ 
bris around M51. At higher surface brightness (pb< 26.5) 
we see many well-known features, including M51’s south¬ 
ern tidal tail, the east and west tidal streams emanating from 
M51b’s disk, and the three-pronged structure (the “crown”) 
just north of M51b. At extremely low surface brightness, we 
trace the extended Northwest Plume, and identify two new, 
extremely faint features—the Northeast and South Plumes. 

We also show in Figure [3] a comparison between our optical 
image and the H I data from Rots et al. (1990). 

The most extended feature is the Northwest Plume. Orig¬ 
inally detected by |Burkhead| (| 1978j), its tip lies 20' (43 kpc) 
from the center ofM51a. While it partially blends with the 
brighter West Streams extending from M51b, a variety of ar¬ 
guments indicate it is a separate feature. First, the Plume 
is morphologically distinct—it extends almost directly west, 
while the northern West Stream extends first northwest and 
then curves toward the northeast (FigurefTli. The Northwest 
Plume is also ~0.06 magnitudes reddeirjthan the adjacent 


West Streams (Table [T] we note that the plume shows a uni¬ 
form color from end to end), impl ying d istinct (albeit similar) 
stellar populations. Durrell et al. ( 2003] ) also note bimodality 
in the PNe velocities northwest of M51 (although not cospa- 
tial with the Plume), arguing for kinematically distinct popu¬ 
lations. Indeed, we note a possible third component, a slight 
protrusion just to the northwest of M51a (marked with an ar¬ 
row in Figure [Tb which may be a stellar stream overlapping 
M51a’s outer disk. The THINGS H I map of M51 ( [Walter et[ 
also shows gas extending roughly in this direction. 


here is very little high column density H I associated with 
the Northwest Plume (Figure [3J, but evidence exists for ex¬ 
tremely diffuse, high-velocity gas in the vicinity (Appleton 
et al.|1986; Rots et al.|1990| Pisano private communication). 


Diffuse H I associated with the Northwest Plume would also 
imply an origin with the gas-rich spiral M51a, as the compan¬ 
ion galaxy M5 lb appears to be extremely gas-poor. 

Dynamically, the Northwest Plume appears several hundred 
Myrs old; its linear extent (~40 kpc in projection) would re¬ 
quire an unrealistic expansion velocity if, for example, it had 


formed during the most recent passage pro posed by [S alo & 
|Laurikainen| ( |2000| (50 to 100 Myr ago) or Howard & Byrd 


2 While this difference is comparable to our absolute photometric uncer- 


tainty, relative color differences between adjacent features are constrained to 
~ 0.01 mag uncertainty, due to the stability of the background corrections. 
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Figure 3. Comparison between our optical image and (in green contours) the H I data ofjRots et al.|jl990). 


(1990) (70 Myr ago). At the rotation speed of M51a (210 
km/s;(Appleton et al.|l986f, a timescale of ~200 Myr is more 
likely. We have also used stellar population modeling (jKotulla 


et al. 20091 to show the Plume is dominated by old stars. Ex¬ 


ponentially declining star formation histories typical of Sb/Sc 
galaxies (r =5-9 Gyr) produce colors too blue (B - V < 0.6) 
to match the Plume; its red color demands much older star 
formation histories (r =3 Gyr). Similarly, any subsequent 
interaction-induced star formation in the Plume must be weak 
(<5% of the Plume’s stellar mass), else the colors would too 
blue by > 0.15 mags even after 500 Myr. As dust extinction 
appears negligible due to the low gas column density (Apple- 
ton et a l.|1986[ |Rots et al.|1990) , our inference of old stellar 
populations appears sound. 

Partially overlapping with the Northwest Plume, but with 
contrasting morphology and color, are the East and West 
Streams. Both streams are bifurcated, and all branches have 
very similar colors (B-V - 0.73), which are in turn simi¬ 
lar to M51b’s average color (B-V = 0.7; Lee et al.|2012[ i. 
Given this, and the streams’ symmetry around and proximity 
to M51b, it seems likely that these streams originated from 
that galaxy’s disk. A variety of processes may give rise to bi- 
furcated tida l streams, includin g kinematic caustics (|S truck & 


Smith [20l2( Smith et ~aI7||2014 > and differences in th e kine- 


matic evolution of the gas and stellar components ( |Mihos| 
[2001) . In the latter mechanism, stellar bifurcations would 
only arise via subsequent star formation in the displaced gas 
tail, which would then show younger stellar populations than 
the collisionless tail. The lack of H I and the uniformly red 


colors of these streams would seem to argue against such a 
scenario in M51b, and instead favor collisionless kinematic 
caustics. 

Based on the colors shown in Figure [2] much of the crown 
north of M51b also appears associated with the compan¬ 
ion, excepting the bluer middle prong. Lee et al.j ( |2012) 
showed that the latter is an extension of M51a’s northern tidal 
arm; indeed, its color is similar to that of the Southeast Tail 
(B-V=0.64), the blueness of which is expected from its ori¬ 
gin in M51a’s disk. We also note the incredible sharpness of 
the ridge formed by the crown’s eastern prong, indicative of a 
strong caustic. 

In the Southeast Tail, the clear separation between the gas 
and stars (Figure [3] noted previously by [Howard & Byrd| 
(|l990j) and (Rots et al.j (Jl990|) is of particular interest. The 
stellar tail curves sharply northward and disappears into the 
tidal detritus near the East Streams, while the more loosely 
wrapped H I tail is much more extended and contains no dif¬ 
fuse stellar counterpart down to our limit of /jb- 30.5. Simu¬ 
lations have shown that such gas/star offsets in tidal tails can 
arise from differences in the radial extent of the two compo¬ 
nents in the pre-encounter disk ( Mihos|2001 1, where the stel¬ 
lar component of a tail consists preferentially of material with 
low specific angular momentum and thus “trails” the high an¬ 
gular momentum gaseous tail. While this may explain the 
tail’s lack of an old stellar population, the lack of either ongo¬ 
ing star formation ( jGil de Paz et al.|2007 ) or a young stellar 
population is more challenging, as the compressive nature of 
tidal interactions is known to form massive star clusters and 
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tidal dwarf galaxies (Barnes & Hcrnquist| 1992; [Jarrett et al. 
2006) in many tidal tails. The lack of a young population in 
M51 ’s tail may indicate a highly inclined encounter, where 
compressive forces in the tail are greatly reduced. As shown 
by Rots et al7[ (jl990J>, an inclined encounter also explains the 
tail’s observed cou nter-rotation, and is re produced well by in¬ 
teraction models of |Salo & Laurikainen| ( |2000| . 

Aside from the structure of the tidal features, our data also 
shed light on the structure of M51a’s outer disk. The colors 
get redder by 0.2-0.4 magnitudes outside of the spiral arm 
region, visible in the colormap of Figure [2] and explicitly in 
Figure |4] which shows photometric profiles along an ellipti¬ 
cal wedge which follows the curvature of the outermost spiral 
arm westward into the outer disk. We see a dramatic drop 
in surface brightness between 220" and 250", corresponding 
precisely with a change in the B-V color from 0.46 to 0.65 
(roughly, the integrated colors of an Sc and Sa galaxy, re¬ 
spectively; |Roberts & Haynes||1994|). This region marks the 
edge of the outermost spiral arm; beyond, the redward gra¬ 
dient continues, ultimately reaching colors similar to those 
of M51b. The dramatic drop in surface brightness is simi¬ 
lar to the arm-interarm contrast measured in the same region 
in M51a by Schweizer ( 1976), and seen in the outer regions 
of other spiral galaxies (ISchweizerll 19761 ITacconi & Young 
[TWO) . 

The red colors, smooth surface brightness profile, and 
round outer isophotes in this region might be interpreted as 
marking the transition to a kinematically hot spheroidal com¬ 
ponent, but we reject this notion for a variety of reasons. 
M51a lacks a strong bulge ((B:T) 7 = 0.27; Pompei & Na- 
tali[|1997) , arguing that the outer component is not likely to 


be extended bulge light. Furthermore, such a bright halo 
component is unlikely given the radius (~ 10 kpc) and sur¬ 
face brightness (p B ~ 25) of the component; by comparison, 
M31’s smooth halo has the much fainter surface brightness of 
p B ~ 2 9-30 at lOkpc (extrapolating from the power law fits 
of Ibata et al. 2014 1 . The most reasonable interpretation for 
this component is that we are simply seeing the outer, redder 
stellar disk in M51a. 

An outer disk origin for this light is also suggested by its 
similarity in color to the more extended Northwest Plume. 
The linear structure of the Plume would require a dynamical 
origin from material with coherent motion and high angular 
momentum, as would be found in an outer disk. Such disks 
are also expected to show redder colors due to effects from 
radial migration (e.g. |Roskar et al.|2008), and such colors are 
frequently observed in outer disks (e.g. Bakos et al. 2008 who 
found typical outer disk colors of g - r ~ 0.5-0.6 = 0.7-0.8 
in B-V). Hence, it seems likely that the red outer regions of 
M51a are simply part of its outer disk, and that the Northwest 
Plume may have been ejected from the outer disk during the 
interaction. 

3.2. Comparison to simulations 

Simulations of the M51 system fall broadly into two 
classes: parabolic single- pass m odels with typical timescales 
of > 300Myr ( Toomre & Toomre|1972| Hernquist|1990 Dur- 

rell et al. 2003) , and more tightly bound, multiple passage 

models with one encounter se veral hundred Myr ago, and a 
second passage <100 Myr ago ( Howard & Byrd 1990| Salo & 


Laurikainen 2000 1 Dobbs et al. 2010). In general, these mod 

els were developed to match the structure and kinematics of 
the H I and the bright tidal debris (the East and West Streams 
and Southern Tail). The degree to which they successfully 


reproduce the fainter, more extended debris (the Northwest, 
Northeast, and South Plumes) is less clear. Here we develop 
new insights gleaned from a comparison between our deep 
imaging data and the existing simulations. 

Perhaps most useful is a comparison involving the North¬ 
west Plume. Similar features arise in simulations of M51 from 
the material nearest the companion at pericenter, which re¬ 
ceives the most dispersive tidal kick and is ejected as a broad 
tidal plume (in contrast, material on the opposite side of the 
disk forms a coherent tidal tail). Tidal features resembling 
the Northwest Plume appear in a variety of M51 simulations, 
but single-passage models tend to reproduce the Plume’s mor¬ 


phology better than the best-fit multipassage models of Salo 
& Laurikainen (2000) or Dobbs et al. (2010): see, for exam- 
ple. Figures 1 and 8 of (Salo & Laurikainen 2000. Howard 

we note that both |Salo 
& Laurikainen ( 2000) and [Dobbs et al.| ( j2010) elaborate on 


only one such model; different parameters in a multi-passage 
encounter may reproduce the feature more effectively. While 
the Northwest Plume cannot unambiguously differentiate be¬ 
tween the model classes, the dynamical age implied by its lin¬ 
ear extent argues that the initial passage of this encounter (or 
the sole passage) must have taken place at least a few hundred 
Myr ago (in more agreement with|Howard & Byrd|1990). 

Unfortunately, the other faint features—the Northeast and 
South Plumes—have no clear analogues in any published sim¬ 
ulation. The Northeast Plume perhaps most resembles the dif- 
fus e northern component seen in the multiple passage models 
by |Salo & Laurikainen ( [2000 i. Single passage models gen¬ 
erally completely fail to reproduce this northern componen 


(Salo & Laurikainen 


Durrell et al. 


component 

__ __, so the pres¬ 
ence of both the crown and this Northeast flume seem to favor 
the multiple-passage models. Indeed, this feature may simply 
be a faint outer extension of material in the crown itself; its 
low luminosity (~ lO 6 L 0r g) would place it below the mass 
resolution of most published simulations, and could explain 
why they miss this feature. 

The South Plume, on the other hand, shows no morpholog¬ 
ical continuity with any adjacent feature, either in optical or 
in H I, and has_ no counterpart in any simulation published to 
date. While Appleton et al. ( j!986 i detected some diffuse H I 
on the galaxy’s southwestern side, it does not overlap with 
the South Plume itself. Given all this, the Plume would seem 
to be either a distinct tidal feature not replicated in previous 
dynamical modeling (although like the Northeast Plume, it 
may be too low mass to discern in published simulations), or 
the remains of a disrupted low luminosity dwarf galaxy in the 
system. 

Regarding the Southeast H I tail, the lack of an old, co- 
spatial stellar component may argue that the gas came from 
an extended gas-rich disk belonging to M51a. Simulations of 
the M51 system favor a highly inclined encounter to explain 
the ki nematics of the tail (Rots et al. 1990, Salo & Laurikainen] 
2000); such a configuration would suppress caustic formation 


in the tail and explain its paucity of star formation and young 
stars. However, the gas/star offset observed in the tail is still 
problematic, as such features should be most prominent in 
highly prograde encounters where the orbital and rotational 
angular momentum is well-aligned ( |Mihos||2001) . Thus the 
data still leave some ambiguity in the dynamical evolution of 
the Southeast Tail. 

In summary, no single model of M51 provides a perfect 
match to the additional constraints derived from our imag¬ 
ing. The single-passage models appear to better reproduce 
















































































































6 


Watkins et al. 



r ( arcsec ) 


Figure 4. Surface brightness and color profiles along the wedge highlighted in the inset, plotted as a function of the average radial distance from M5 la’s center. 
Error bars in both plots are calculated using the background variance, and are often smaller than the point size. 


the diffuse Northwest Plume, while the multi-passage models 
reproduce better the diffuse northern material (the crown and 
Northeast Plume). No simulations to date have produced any¬ 
thing resembling the South Plume, although its low luminos¬ 
ity suggests that higher resolution simulations may be neces¬ 
sary to model its structure. The lack of a stellar counterpart in 
M51’s Southeast HI tail argues that the tail was drawn almost 
exclusively from the outermost, gas-rich regions of M51a’s 
disk, or that extensive kinematic segregation of stars and gas 
has occurred, without triggering any wide-spread star forma¬ 
tion in the tidally stripped gas. 

Future modeling campaigns thus are warranted, and have a 
variety of new constraints and challenges ahead. New simu¬ 
lations must recover the newly discovered tidal plumes, en¬ 
sure that the Northwest Plume is correctly oriented and gas- 
free, reproduce the strong offset between the stars and gas 
in the Southeast Tail, and bifurcate M51b’s East and West 
Streams. Successful simulations will then be able to pin down 
the timescale and geometry of the encounter—critical inputs 
to studies of spiral structure, ISM evolution, and star forma¬ 
tion in this classic interacting system. 


The work was supported by the NSF through grant 1108964 
to JCM. 
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